Bianchi type I cosmology with scalar and spinor fields 



Bijan Saha and Todor Boyadjiev 

Laboratory of Information Technologies 
Joint Institute for Nuclear Research, Dubna 
141980 Dubna, Moscow region, Russi^ 

We consider a system of interacting spinor and scalar fields in a gravitational 
field given by a Bianchi type-I cosmological model filled with perfect fluid. The in- 
teracting term in the Lagrangian is chosen in the form of derivative coupling, i.e., 
=^nt = J ^>,a^>'^F 1 with F being a function of the invariants I an J constructed from 
bilinear spinor forms S and P. We consider the cases when F is the power or trigono- 
metric functions of its arguments. Self-consistent solutions to the spinor, scalar and 
BI gravitational field equations are obtained. The problems of initial singularity and 
asymptotically isotropization process of the initially anisotropic space-time are stud- 
ied. It is also shown that the introduction of the Cosmological constant (A-term) 
in the Lagrangian generates oscillations of the BI model, which is not the case in 
absence of A term. Unlike the case when spinor field nonlinearity is induced by 
self-action, in the case in question, wehere nonlinearity is induced by the scalar field, 
there exist regular solutions even without broken dominant energy condition. 
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1. INTRODUCTION 

Nonlinear generalization of classical field theory remains one of the possible ways to 
overcome the difficulties of the theory which considers elementary particles as mathemat- 
ical points. The gravitational field equation is nonlinear by nature and the field itself 
is universal and unscreenable. These properties lead to definite physical interest for the 
proper gravitational field to be considered. Nonlinear self-couplings of the spinor fields 
may arise as a consequence of the geometrical structure of the space-time and, more pre- 
cisely, because of the existence of torsion. As early as 1938, Ivanenko [l|, 0, 0| showed 
that a relativistic theory imposes in some cases a fourth-order self-coupling. In 1950, 
Weyl |3] proved that, if the affine and the metric properties of the space-time are taken 
as independent, the spinor field obeys either a linear equation in a space with torsion or 
a nonlinear one in a Riemannian space. As the self-action is of spin-spin type, it allows 
the assignment of a dynamical role to the spin and offers a clue about the origin of the 
nonlinearities. A nonlinear spinor field, suggested by the symmetric couplingbetween 
nucleons, muons, and leptons, has been investigated by Finkelstein et. al. p in the 
classical approximation. Although the existence of spin- 1/2 fermion is both theoretically 
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and experimentally undisputed, these are described by quantum spinor fields. Possible 
justifications for the existence of classical spinors has been addressed in 

The present cosmology is based largely on Friedmann solutions of the Einstein equa- 
tions, which describe the completely uniform and isotropic universe ( "closed" and "open" 
models) The main feature of these solutions is their non-stationarity. The idea of an 
expanding Universe, following from this property, is confirmed by the astronomical ob- 
servations and it is now safe to assume that the isotropic model provides, in its general 
features, an adequate description of the present state of the Universe. Although the Uni- 
verse seems homogeneous and isotropic at present, the large scale matter distribution 
in the observable universe, largely manifested in the form of discrete structures, does 
not exhibit homogeneity of a higher order. Recent space investigations detect anisotropy 
in the cosmic microwave background. In fact, the theoretical arguments j3| and recent 
experimental data that support the existence of an anisotropic phase that approaches 
an isotropic one leads to consider the models of universe with anisotropic back-ground. 
Zel'dovich was first to assume that the early isotropization of cosmological expanding pro- 
cess can take place as a result of quantum effect of particle creation near singularity . 
This assumption was further justified by several authors [1, III- 

A Bianchi type-I (BI) Universe, being the straightforward generalization of the flat 
Friedmann-Robertson- Walker (FRW) Universe, is of particular interest because it is one 
of the simplest models of an anisotropic Universe that describes a homogeneous and 
spatially flat Universe. Unlike the FRW Universe which has the same scale factor for each 
of the three spatial directions, a BI Universe has a different scale factor in each direction, 
thereby introducing an anisotropy to the system. It moreover has the agreeable property 
that near the singularity it behaves like a Kasner Universe, even in the presence of matter, 
and consequently falls within the general analysis of the singularity given by Belinskii et 
al Also in a Universe filled with matter for p — C,£, ^<l,it has been shown that 
any initial anisotropy in a BI Universe quickly dies away and a BI Universe eventually 
evolves into a FRW Universe Since the present-day Universe is surprisingly isotropic, 
this feature of the BI Universe makes it a prime candidate for studying the possible effects 
of an anisotropy in the early Universe on present-day observations. 

It should be noted an important property of the isotropic model is the presence of a 
singular point to time in its space-time metric which means that the time is restricted. 
Is the presence of a singular point an inherent property of the relativistic cosmological 
models or it is only a consequence of specific simplifying assumptions underlying these 
models? Motivated to answer this question self-consistent svstem of nonlinear spinor and 
BI gravitational fields were studied in a series of papers [l^, EE EE E3- should be 
mentioned that a spinor field in a BI universe was also studied by Belinskii and Khalat- 
nikov Using Hamiltonian techniques, Henneaux studied class-A Bianchi universes 

generated by a spinor source 

In this report we consider a self-consistent system of spinor, scalar and Bianchi type-I 
gravitation fields in presence of perfect fluid and cosmological constant A. It should be 
noted that the inclusion of the A adds new dimension in the evolution of the universe. 
Assuming that the A term may be both possitive and negative, it opens a much wider 
field of possibiliteis in the search for a singularity-free solutions to the field equations. 
Apart from the papers written on this subject by us earlier, in this one beside the power 
type interaction term we consider a trigonometric one as well. In addition to asymptotic 
analysis, numerical analysis of the corresponding nonlinear differential equations has been 
performed to some extent as well. 
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2. DERIVATION OF BASIC EQUATIONS 

In this section we derive the fundamental equations for the interacting spinor, scalar 
and gravitational fields from the action and write their solutions in term of the voulume 
scale T defined bellow fl2.27|) . We also derive the equation for T which plays the central 
role here. 

We consider a system of nonlinear spinor, scalar and BI gravitational field in presence 
of perfect fluid given by the action 

y{g; Y, \jf,(p) = l ^V^dQ. (2.1) 

with 

The gravitational part of the Lagrangian ()2.2|) is given by a Bianchi type I (BI hereafter) 
space-time, whereas the remaining parts are the usual spinor and scalar field Lagrangian 
with an interaction between them and a perfect fluid as well. 



A. Material field Lagrangian 

For a spinor field i//, symmetry between \}/ and \ff appears to demand that one should 
choose the symmetrized Lagrangian jl^l. Keep it in mind we choose the spinor field 
Lagrangian as 

with m being the spinor mass. 

The mass-less scalar field Lagrangian is chosen to be 



m\]/\]/, (2.3) 



^sc = ^(p,a<. (2.4) 
The interaction between the spinor and scalar fields is given by the Lagrangian [1 

^ni=^(P,a<F, (2.5) 

with A being the coupling constant and F is some arbitrary functions of invariants gen- 
erated from the real bilinear forms of a spinor field with the form 

S = iffXi/ (scalar), (2.6a) 

P = i^ff^r (pseudoscalar), (2.6b) 

= {xj/j'^Xj/) (vector), (2.6c) 

= (xfry^y^xir) (pseudovector), (2.6d) 

Q^"^ = (^a'^^v^) (antisymmetric tensor), (2.6e) 
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where <7^^ = {i/2)[y^'Y^ — y^y^]. Invariants, corresponding to the bihnear forms, are 



/ = 5^ (2.7a) 

J = P^, (2.7b) 

Iv = v^v^ = {xlry^xir)g^,y{xlry''xif), (2.7c) 

Ia = A^A^ = {ifry^f'v)8nv{¥ffv), (2.7d) 

Iq = QiivQ^'' = (V/C7'^»g^a^v/3(V^C7«V)- (2.7e) 



According to the Pauh-Fierz theorem [22| among the five invariants only / and J are 
independent as all other can be expressed by them: ly = —Ia = I + J and Iq = I — J. 
Therefore, we choose F —F{I,J), thus claiming that it describes the nonlinearity in the 
most general of its form Note that setting A = is (12. 5|) we come to the case with 
minimal coupling. 

The term describes the Lagrangian density of perfect fiuid. 



B. The gravitational field 

As a gravitational field we consider the Bianchi type I (BI) cosmological model. It is 
the simplest model of anisotropic universe that describes a homogeneous and spatially fiat 
space-time and if filled with perfect fiuid with the equation of state p = ^e, C < 1; it 
eventually evolves into a FRW universe The isotropy of present-day universe makes 
BI model a prime candidate for studying the possible effects of an anisotropy in the early 
universe on modern-day data observations. In view of what has been mentioned above 
we choose the gravitational part of the Lagrangian ()2.2|) in the form 

i?, = ^. (2.8) 

where R is the scalar curvature, fc = SnG being the Einstein's gravitational constant. The 
gravitational field in our case is given by a Bianchi type I (BI) metric 

ds^ = df-a^dx^-b^dy^-c^dz^, (2.9) 

with a, c being the functions of time t only. Here the speed of light is taken to be unity. 
The metric ()2.9|) has the following non-trivial Christoffel symbols 



Tio - -, r^o--, r3o-- 



r?i = ad, = bh, = cc. 



(2.10) 
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The nontrivial components of the Ricci tensors are 



R 



'a b c 
-(- + - + - 
.a b c 

d d fb c 
- + -- + - 
a a\b c 



R 



R 





\b b 


fc d 






[- + - 






\c a 






^d b 










Lc c ^ 


\a b 



From ()2.11|) one finds the following Ricci scalar for the BI universe 

^/d b c db be cd 
i?=-2 - + - + - + -- + -- + — 
\a b c ab be c a 

The non-trivial components of Riemann tensors in this case read 



R 



,01 
01 



^ p02 



b 
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c 

1 

c 



R 



12 
12 



(2.11a) 
(2.11b) 
(2.11c) 
(2.11d) 

(2.12) 



(2.13) 



,.23 _ be p3i _ c« 
ab be ea 

Now having all the non-trivial components of Ricci and Riemann tensors, one can easily 
write the invariants of gravitational field which we need to study the space-time singularity. 
We return to this study at the end of this section. 



C. Field equations 

Let us now write the field equations corresponding to the action ()2.ip . 
Variation of 1)2.111 with respect to spinor field Yiw) gives spinor field equations 

iy^'^uMf-mxif + ^x^f + ^iy^x^/ = 0, (2.14a) 
iVnWr^ +m\ff-^\j/-':^i\j/"^ = 0, (2.14b) 

where we denote 

Since the nonlinearity in the spinor field equations is generated by the interacting scalar 
field, the Eqs. ()2.14|) can be viewed as spinor field equations with induced nonlinearity. 
Varying ()2.H) with respect to scalar field yields the following scalar field equation 



1 



■(v^^"'^(1+AF)9,m)=0- (2-15) 



Finally, variation of ()2.1|) with respect to metric tensor g^y gives the Einstein's field 
equation in account of the A-term has the form 

Rl-\d; = KT;-5;A. (2.16) 
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ktI-a, 


(2.17a) 


KTi-A, 


(2.17b) 


kT,'-A, 


(2.17c) 


kT^-A, 


(2.17d) 



In view of 1)2.111) and ()2.12|) for the BI space-time ()2.9|) we rewrite the Eq. ()2.16|) as 

b c be 

b c be 

c d c d 

- + -H 

c a c a 

d b db 

a b ab 
db be c d 
ab be c a 

where over dot means differentiation with respect to t and Ty is the energy-momentum 
tensor of the material field given by 

Here T^^^ is the energy-momentum tensor of the spinor field defined by 

T4i,= '-^gp''(^yi^yvW+wyvypiW-^iiWyvW-^vwyiiW^ -s^^sp. (2.19) 

The term with respect to ()2.14j) takes the form 

^sp = -[&S + '^P). (2.20) 
The energy-momentum tensor of the scalar field T^^ii is given by 

Tj^ = W^-5;^,,. (2.21) 

For the interaction field we find 

TZ^ = lFcp^^cp^^-5l^^,. (2.22) 

^^{m) energy-momentum tensor of a perfect fluid. For a universe filled with perfect 

fluid, in a comoving system of reference such that m'^ = (1, 0, 0, 0) we have 

T^^im) = {P + £>nU^ - 5> = (e, -P^ -P^ -P)^ (2.23) 

where energy of the perfect fluid £ is related to its' pressure p by the equation of state 
p = C,£. Here C, varies between the interval < ^ < 1, whereas ^ = describes the 
dust Universe, C = ^ presents radiation Universe, ^ < ^ < 1 ascribes hard Universe and 
C, = \ corresponds to the stiff matter. 

In the Eqs. fl2.14|l and ()2.19|) is the covariant derivatives acting on a spinor field as 

mm 

^MV^ = 0-rMV^, v^v^=0 + v^r^, (2.24) 

where are the Fock-Ivanenko spinor connection coefficients defined by 

^li = \f(^layv-di,ya). (2.25) 
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For the metric ()2.9p one has the following components of the spinor connection coefficients 
ro = 0, T, = \d{t)tf, T2 = \b{t)ff, T, = \c{t)ff. (2.26) 

The Dirac matrices y^{x) of curved space-time are connected with those of Minkowski 
one as follows: 

f = f, i = fla, f = f/b, ^ = f/c 



with 



f-i'^-X -^-f-r-ro 



where a,- are the Pauli matrices: 

^ =^10j' ^ J' ^ =^0-1 

Note that the y and the O matrices obey the following properties: 

fyj + yjf = 2r]'J^ /,j = 0,l,2,3 
fy5 + 'p5f = 0, {ff=I, / = 0,1,2,3 
o^o^ = 5jk + iejMo\ J, ^, / = 1 , 2, 3 

where rj/j = {1, — 1, — 1, — 1} is the diagonal matrix, dj^ is the Kronekar symbol and Ejki is 
the totally antisymmetric matrix with £123 = +1. 

We study the space-independent solutions to the spinor and scalar field equations ()2.14j] 
and ()2.15|) so that 1//= and (p = (p{t). Defining 

i^abc^ (2.27) 

from ()2.15|) for the scalar field we have 

The spinor field equation (j2.14aj) in account of (j2.24j) and (j2.26j) takes the form 



-■KM) 



VZ-rnvz+^VZ + ^/y^V^-O. (2.29) 



Setting Vj{t) = y/t\ffj{t), j = 1,2,3,4, from ()2.29|1 one deduces the following system of 
equations: 

Vi + i{m-^)Vi-m3 = 0, (2.30a) 
V2 + iim-^)V2-':^V4 = 0, (2.30b) 
V3~i{m-^)V3 + mi = 0, (2.30c) 
V4-i{m-^)V4 + ':^V2 = 0. (2.30d) 
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From ()2.14a|) we also write the equations for the invariants 5, P and A = i/zy^)^!// 

So-2^Ao = 0, (2.31a) 
Po-2(m-^)Ao = 0, (2.31b) 
Ao + 2(m-^)/'o + 2^5o = 0, (2.31c) 

where 5o = tS, Pq = tP, and Aq = xA. The Eqn. fl2.31|) leads to the following relation 

S^+P^+A^ = C^/t^, = const. (2.32) 

Giving the concrete form of F from ()2.3()|1 one writes the components of the spinor 
functions in explicitly and using the solutions obtained one can write the components of 
spinor current: 

/ = xj/yt^y/. (2.33) 

The component 

f = - \v*Vi + ViV2 + V;V3 + VM , (2.34) 

T 

defines the charge density of spinor field that has the following chronometric-invariant 
form 

P = (J0-/)^/'. (2.35) 
The total charge of spinor field is defined as 

oo 

Q= J py^^dxdydz = ptr, (2.36) 

— oo 

where is the volume. From the spin tensor 

5^^'^ = iv^{f a^Va^Y}v^- (2-37) 
one finds chronometric invariant spin tensor 

^cf =(%o5'^'")'^'' (2-38) 
and the projection of the spin vector on k axis 

oo 

Sk= J S'J^f^/^dxdydz = S'J^\v. (2.39) 

— oo 

Let us now solve the Einstein equations. To do it we first write the expressions for the 
components of the energy-momentum tensor explicitly: 



(2.40) 



2t2(i + AF) 
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In account of ()2.40|) subtracting ()2.17aj) from ()2.17b|) . one finds the following relation 
between a and b 

" - (2.41) 



Analogically, one finds 



a ^ f dt 



b ^ fdt 
-=D,e.viX,J- 



[2A2) 



Here Di, D2, D^, Xi, X2,X2 are integration constants, obeying 

DiD2D3 = l, Xi+X2+X3 = 0. (2.43) 
In view of ()2.43|) from ()2.41|) and ()2.42j) we write the metric functions explicitly ^ 



a(t) 
b{t) 
c{t) 



(D?D3)i/3Ti/3exp 



2X1 +X3 



dt 



3 J x{t) 
X1-X3 f dt 



{D,Dl 



-1/3^1/3 



exp 



3 J x{t 
X1+2X3 r dt 



T(0 



(2.44a) 
(2.44b) 
(2.44c) 



As one sees from (I2.44a|l . (j2.44bjl and ()2.44c|l . for T = with n > 1 the exponent tends to 
unity at large t, and the anisotropic model becomes isotropic one. 

Further we will investigate the existence of singularity (singular point) of the grav- 
itational case, which can be done by investigating the invariant characteristics of the 
space-time. In general relativity these invariants are composed from the curvature tensor 
and the metric one. In a 4D Riemann space-time there are 14 independent invariants. 
Instead of analyzing all 14 invariants, one can confine this study only in 3, namely the 
scalar curvature I\ =R, h = ^^vM^' ^^"^ Kretschmann scalar 73 =Raj5^v^'^^^^ ■ 
regular space-time point, these three invariants Ii, 12,13 should be finite. Let us rewrite 
these invariants in detail. 

For the Bianchi I metric one finds the scalar curvature 



^/d b c db be cd 
\a b c ab be c a 



Since the Ricci tensor for the BI metric is diagonal, the invariant I2 
a sum of squares of diagonal components of Ricci tensor, i.e.. 



(2.45) 

R^yR^''=RlR^ is 



h 



{R'^'+{R\Y+{R'2)'+{Rl) 



,U2 



2\2 



,3X2- 



(2.46) 



with the components of the Ricci tensor being given by (j2.1ip . 

Analogically, for the Kretschmann scalar in this case we have 73 
of squared components of all nontrivial 7?^^V) which in view of ()2.13p can be written as 

h = 



R 



01 



+ (^2 



+ (R'\3 



= 4 



a\2 /b\^ /c\^ /db\'2 /bc\'2 /cd\^ 
a/ \b/ \c/ \ab/ \bc/ \ca/ . 



(2.47) 
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Let us now express the forgoing invariants in terms of T. From Eqs. ()2.44|) we have 



ai = A,T^/3exp^(y,/3) Jt^^d^j, (2.48a) 
(/= 1,2,3,), (2.48b) 



ai 3 T 

a, ^ (y, + i)(y,-2) 1 

ai 9 t2' 



(2.48c) 



i.e., the metric functions a,b,c and their derivatives are in functional dependence with T. 
From Eqs. ()2.48|) one can easily verify that 

1 1 1 

Thus we see that at any space-time point, where T = the invariants h, 12,13, as well as 
the scalar and spinor fields become infinity, hence the space-time becomes singular at this 
point. 

In what follows, we write the equation for t and study it in details. 
Summation of Einstein equations ()2.17a|) , ()2.17b|) , ()2.17cp and ()2.17d|) multiplied by 3 
gives 

I = ^K(^m5 + ^5 + 5fP + e-p) -3A. (2.49) 

For the right-hand-side of ()2.49|) to be a function of T only, the solution to this equation 
is well-known psl ]. 

Let us demand the energy- momentum to be conserved, i.e., 

T;^;y = 7;r,v +r^v7? -rJvTp' = 0, (2.50) 
which in our case has the form 

'-{rTSy-'-Tt-tr^-iTi^O. (2.51) 

In account of the equation of state p = and 

{m-&)So-^Po = 
which follows from ()2.31|) . after a little manipulation from ()2.51|) we obtain 

e = eo/T'+^, p=^£o/t'+^. (2.52) 
In view of ()2.52|) the Eq. ()2.49|) can be written as 

^ = ^K(mS+^S + ^P+{l - C)eo/T^+^) -3A. (2.53) 

As it was mentioned earlier, we consider F as a function of I, J or I ± J. In the section to 
follow we study the Eq. (j2.53|) in details. 
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3. EXACT SOLUTIONS AND NUMERICAL ANALYSIS 

In the preceding section we solved the spinor, scalar and gravitational field equations 
and wrote the solutions in terms of volume-scale T. It was also mentioned that for the 
right hand side of the Eq. ()2.53|) to be the function of T, this equation is quadratically 
integrable. In what follows, we explicitly write the solutions for corresponding equations 
given some concrete form of T. 



A. Exact solutions 

Here we consider the cases with minimal coupling and with F being the function of 
either I or J (with zero mass). In this subsection we simply write the solutions to the 
spinor field equations explicitly and present the solution for T in quadrature. 



1. Minimally coupled scalar and spinor fields 

Let us first consider the case with minimal coupling when the scalar and the spinor 
fields interact through gravitational one. In this case from (j2.31|) one finds S = Cq/x. 
Scalar field and the components of the spinor field in this case have the following explicit 
form 

(p = c/^, (3.1) 



V ^ 



with the integration constants Cj satisfying Cq = Cj -|-C| — C| — C|. 
Eq. ()2.53p in this case takes the form 



with the solution 



(3.2) 



x=-K(,nCo + eo{\-C)/T^)-3Ax, (3.3) 



V3t. (3.4) 



K(mCoX + £oX^-^^ -Ax^ + E 

Here E is the constant of integration. Let us note that being the volume-scale x cannot 
be negative. On the other hand the radical in ()3.4|1 should be positive. This fact leads to 
the conclusion that for a positive A the value of x is bound from above giving rise to an 
oscillatory mode of expansion of the BI space-time. 
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2. Case with F = F(/) 

Here we consider the interacting system of scalar and spinor field with the interaction 
given by = (A/2)(p^9^F(/). As in the case with minimal coupling from ()2.31aj) one 
finds 

5 = —, Co = const. (3.5) 

T 

For components of spinor field we find [T 



(3.6) 



with Ci being the integration constants and are related to Cq as Cq = C\ ^ C\ — C\ — C\. 
Here [5 = J{m-^)dt. 

For the components of the spin current from ()2.33j) we find 

/ = L[cj+Cj + Cj + Cl], / = ^[CiC4 + C2C3]cos(2/3), 
/ = ^[CiC4-C2C3]sin(2/3), / = A [QCg -C2C4]cos(2i8), 
whereas, for the projection of spin vectors on the X, Y and Z axis we find 
j-23,0 ^ C1C2 +C3C4 ^31,0 ^12,0 ^ — C| + C| — C4 

bet ' ' 2abx 

Total charge of the system in a volume 7^ in this case is 

Q=[C^i+Cj + Cf + C|] y. (3.7) 

Thus, for T 7^ the components of spin current and the projection of spin vectors are 
singularity-free and the total charge of the system in a finite volume is always finite. 
The equation for determining T in this case has the form 

f=-K(mCo + n + eo{l-Q/x^^-3AT. (3.8) 

Recalling that ^ = XCoC^Fi/t:^{1 + ?iF{I)f the solution to Eq. (jSSl) can be written in 
quadrature 

/ , = V3t, (3.9) 

fc(^mCoT + C2/2(l + AF) + eoTi-^) -Az^ + E 



with E being some integration constant. Given the explicit form of F{I) for different A 
we have find different mode of expansion. We study this case in details numerically in 
the subsection to follow. 
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3. Case with F = F{J) 

Finally we consider the interacting system of scalar and spinor field with the interaction 
given by = {X /2)(pp^(p^F{J). In the case considered we assume the spinor field to be 
massless. It gives ^ = 0. Note that, in the unified nonlinear spinor theory of Heisenberg, 
the massive term remains absent, and according to Heisenber g, t he particle mass should 
be obtained as a result of quantization of spinor prematter [28|. In the nonlinear gen- 
eralization of classical field equations, the massive term does not possess the significance 
that it possesses in the linear one, as it by no means defines total energy (or mass) of 
the nonlinear field system. Thus without losing the generality we can consider mass-less 
spinor field putting m — Q. Then from (j2.31bjl one gets 

P = Dq/x, Do = const. (3.10) 

In this case the spinor field components take the form 



(3.11) 



The integration constants D,- are connected to Dq by Dq = 2 (Z)j -\-D\ — D3 — D|). Here we 
set a = / '^dt. 

For the components of the spin current from ()2.33p we find 

/ = l[D\+Dl + Dl + D^^, / = — [D2D3+DiD4]cos(2a), 
X ax 

f = ^ [D2D3 -DiD4]sin(2a), / = ^ [D1D3 -D2Z)4]cos(2c7), 

whereas, for the projection of spin vectors on the X, Y and Z axis we find 

bcx ' ' 2abx 

For X in this case we have 

3 



X = ^k{^Co + eo(l - 0/^9 - 3At. (3.12) 



In view of fl3.1()j) . $f in this case takes the form analogical to that taken by ^ in previous 
case with F/ replaced by Fj. Then solution to Eq. ()3.12|) we write in quadrature as 

V3t, (3.13) 



K(^C2/2(l+AF)+eoTi-?) -Ax^+E 



Depending on the form of F and A we have different mode of expansion of BI universe as 
in previous case. In what follows we numerically study the aforementioned cases. 
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B. Numerical experiments 

In this subsection we numerically solve the Eq. (j3.8|) for some different choice of F . As 
it was mentioned earlier, setting A = in ()3.8p we come to the case with minimal coupling 
given by ()3.3|) . whereas, assuming m = we get ()3.12|) . Let us first rewrite Eq. ()3.8p 

T = ^{T,p), (3.14) 

where we denote 

^ = ^K(mCo + ^Co + eo(l-C)/'!^^) -3At, (3.15) 

and p = {k, A,m,Co,C,£o, ^,A} is the set of the parameters. Since in the examples we 
consider F = F{S), let us rewrite ^ in terms of S. In account of 5 = Cq/t for ^ we have 

^ = XC^Fs/2T^{l+XF{S)f. 

From mechanical point of view the Eqn. ()3.14|) can be interpreted as an equation of 
motion of a single par ticle with unit mass under the force ^{T,p). Then the following 
first integral exists j29| 

t = ^/2[E-'^{T,p)]. (3.16) 
Here E is the integration constant and 

3 



K(mCoT + C^/2{l+XF)+£oT-^^ -At2 



is the potential of the force We note that the radical expression must be non-negative. 
The zeroes of this expression, which depend on all the problem parameters p define the 
boundaries of the possible rates of changes of T{t). In what follows we numerically analyze 
fl3.14|l and ()3.15j) for different choice of F(/) as well as for different problem parameters 
P- 



1. F = S" 

Let us first choose F to be a power function of 5(or /), setting F = S'\ In this case 
setting Co = 1 and C = 1 we rewrite ^ as 



^ 3k/ Xnv 

^ = — m + 



with the potential 



2(A + T")2 



(1-C) 



-3At, 



(3.17) 



mx- 



2(A + T«) 



+ eoT 



(3.18) 



Note that the nonnegativity of the radical in ()3.1(i|l in view of ()3.18j) imposes restriction 
on T from above in case of A > 0. It means that in case of A > the value of T runs 
between and some tmaxi where Tmax is the maximum value of T for the given value of p. 
This equation has been studied for different values of parameters p. Here we demonstrate 
the evolution of T for different choice of Tq for fixed "energy" E and vise versa. 
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As the first example we consider massive spinor field with m = 1 . Other parameters 
are chosen in the following way: coupling constant A =0.1, power of nonlinearity n = 4, 
and cosmological constant A = 1/3. We also choose ^ = 0.5 describing a hard universe. 




FIG. 1: Perspective view of the potential 
^ (t) with Bl space-time being filled per- 
fect fluid describing hard universe 



FIG. 2: Evolution of the BI space-time cor- 
responding to the potential given in Fig. ^ 
for different choice of E. 



In Fig. [T]we plot corresponding potential ^ (t) multiplied by the factor 2/3. As is seen 
from Fig. ^ and Fig. |21 choosing the integration constant E we may obtain two different 
types of solutions. For E > 0.5 solutions are non-periodic, whereas for £min 

< £ < 0.5 the 

evolution of the universe is oscillatory. 

As a second example we consider the massless spinor field. Other parameters of the 
problem are left unaltered with the exception of ^. Here we choose C — ^ describing stiff 
matter. It should be noted that this particular choice of ^ gives rise to a local maximum. 
This results in two types of solutions for a single choice of E. 



As one sees from Fig. |21 for E > M there exists only non-periodic solutions, whereas, 
for Emin <E < —0.5 the solutions are always oscillatory. For E G (—0.5, M) there exits two 
types of solutions depending on the choice of Tq- In Fig- IHwe plot the evolution of T for 
E e (— 0.5,M'). As is seen, for Tq G (0, A) we have periodical solution, but due to the fact 
that T is non-negative, the physical solutions happen to be semi-periodic. For Xq G (5, C) 
we again have oscillatory mode of the evolution of T. This two region is separated by a 
no-solution zone {A,B). 

Let us also consider the case with A < 0. For a negative A, as well as in absence of 
the A-term the evolution of T is always exponential as it is seen in Fig. ??. In this case 
the initial anisotropy of the BI space-time quickly dies away and the universe becomes 
isotropic one. 



Let us analyze the dominant energy condition in the Hawking-Penrose theorem j2 



16 



Bijan Saha and Todor Boyadjiev 




FIG. 3: Perspective view of the potential 
(t) with BI space-time being filled with 
stiff matter 



FIG. 4: Evolution of the BI space-time 
corresponding to the potential given in 
Fig. 131 in case of massless spinor field for 
different choice of To with E G (— 0.5,M). 




FIG. 5: Perspective view of T for a negative A. As one 
sees, evolution of the universe in this case takes exponen- 
tial character and the initial anisotropy of the BI space- 
time quickly dies away. 



For a BI universe the dominant energy condition can be written in the form [15[ 

Tq > T^a^ + Tib^ + T^c^, (3.19a) 

To > Tla^, (3.19b) 

To > Tib\ (3.19c) 

To > T^c^- (3.19d) 

It was shown that for the spinor field with self-coupling the regular solutions can be 
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obtained only by violating energy- dominanat condition |15|. Regular solution that does 
not violate this condition was found only for linear spinor case by means of a positive 
A term. To analyze this condition for the spinor field with induced nonlinearity let us 
rewrite the components of energy-momentum tensor. For energy density in this case we 
have 

As one sees from ()3.20j) for any positive value of T energy density is always positive definite. 
As T ^ 0, Tq — > oo, whereas Tq decreases as T increases. For the pressure components we 
have 

= = = Kr^lTcsr I^'^""'" - ') - - (3.21) 

The second term in (j3.21|) is always positive, it means has a greater value when Bl 
universe is filled with dust, i.e., when ^ =0. To investigate the dominant energy condition 
we study the pressure term (since = T2 = T^, hereafter we mention it as T^) in details. 
For simplicity we set C = 1 and Co = 1. It is clear from 1)3.211) that if 

x">X{n-l), (3.22) 

we have < 0. In this case the dominant energy condition remains unbroken. From 
()3.22|) we see, for A = the foregoing inequality holds for any T > 0. It means like the 
linear spinor field, the system of minimally coupled scalar and spinor fields possesses 
regular solutions without violating the dominant energy condition. For an interacting 
system this condition holds for any negative n with a positive A and vice versa. Let us 
now see what happens when both n and A are positive (negative). Note that the coupling 
constant A may take any value. The magnitude of A defines the strength of interaction. 

It is clear that for a large value of T the inequality ()3.22|1 is likely to be held for any 
reasonable value of A. For T close to zero it is unlikely to be held, but in that case, as 
it was mentioned earlier, Tq tends to 00, so the dominant energy condition at this region 
remains safe. Finally we see what happens to it when T is close to unity. Tq at this point is 
reasonably small, whereas for a relatively large n we have a situation when dominates 
Tq. Relative behavior of Tq and T^ has been shown in Fig. IHl Thus we conclude that 
in case of interacting spinor and scalar fields it is possible to construct regular solutions 
without violating dominant energy condition of Hawking-Penrose theorem. 



2. F = sin S 



Let us now consider the case with F = sin S. In this case for ^ we have 
^ 3k/ Acos5 (1-C)\ 



2 V 2T2(A+sin5)2 ^ 



(3.23) 



with the potential 



2(1+ A sin 5) 



+ eoT 



■Ax' 



(3.24) 
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FIG. 6: Perspective view of and Tj for 
a positive n. As one sees for a small value 
of n it is possible to construct a regular so- 
lution without violating dominant energy 
condition. 



FIG. 7: For a large n there exists some 
value of T where the pressure component 
prevails energy. In this case the dominant 
energy condition breaks down. 



It should be noted that unlike the case with F being a power function of T where the 
nonlinearity appears in the region with large value of T, in the case under consideration, 
a number of interesting properties emerge in the region where < T < 1 , namely, in the 
vicinity of the singular point T = 0. A perspective view of the potential ^(t) is given in 
the Figs. IHlandini Here we choose the problem parameters as follows: K = 2/3, spinor mass 
m = 1, coupling constant A = 0.01, cosmological constant A = 2/3, Co = 1 and ^ = 2/3. 




FIG. 8: Perspective view of the potential 
^ (t) with BI space-time being filled per- 
fect fluid describing hard universe 



FIG. 9: Extraordinary behavior of the potential 
in the vicinity of the singular point T = that 
occurs due to the nonlinear term F. 



It is clear from Fig.|Hlthat an oscillatory mode of evolution takes place, as was expected 
for a positive A. 
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Let us now study the system for a negative A. Contrary to the case with F = 5", 
where all the solutions for a negative A grow exponentially, in the case considered here 
an interesting situation occurs for some special choice of parameters. 




FIG. 10: Perspective view of the potential ^{z) with a negative A. 



As one sees from Fig. ??, depending on the integration constant and initial value of T, 
the mode of evolution can be both finite and exponential. For the integration constant 
being at the level AB in Fig. ?? (here it is —3), with Tq G (0, Ta) the evolution of T finite 
and similar to one in Fig. |2l corresponding to £ = 1 , whereas, for Tq > Tg we have T that 
is expanding exponentially. In case of E being at the same level with point M we have 
the similar picture of evolution, but in this case once T reaches point Tm, the process of 
evolution would come to a halt. Thus we conclude that for a trigonometric interaction 
term the system even with a negative A admits non-exponential mode of evolution. 

To investigate the dominant energy condition let us write the components of energy- 
momentum tensor. For simplicity we set Co = 1 and in term of S for energy density we 
write 

Since T is a positive quantity, S is positive as well. As one sees from ()3.25p for any positive 
value of S and A < 1 energy density is always positive definite and proportional to S^. 
Since 5 = 1/t, it means that Tq has its maximum as T — and tends to zero as T — >■ oo. 
On the For the pressure components we have 
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As one sees, for a A < 1 may both positive or negative depending on the sign of 
cos 5. Moreover, its maximum value is proportional to 5^. Thus, in case of F = sin 5 
for all possible values of C, and A (necessarily nontrivial) there exists intervals (5„,5„-|-i) 
such that for S E (SmSn+i) the inequality Tq < takes place as it is shown in Fig. ??. 
Therefore we conclude that the regular solutions obtained in this case results in broken 
dominant energy condition. 




4. CONCLUSION 

Within the framework of the simplest model of interacting spinor and scalar fields it 
is shown that the A term plays very important role in BI cosmology. In particular, it 
invokes oscillations in the model which is not the case when A term remains absent. For a 
non-positive A we find an universe expanding exponentially, hence the initial anisotropy 
of the model quickly dies away, whereas for a positive A with the corresponding choice 
of integration constant S' one finds the oscillatory mode of expansion of the universe. 
In this case it is possible to construct solutions those are always regular. It should be 
emphasized that if the spinor field nonlinearity is generated by self-action the regularity 
of the solutions obtained results in the violation of the dominant energy condition of 
Penrose-Hawking theorem whereas in the case considered here, when the spinor field 
nonlinearity is induced by the scalar one, regular solutions can be obtained even without 
breaking the aforementioned condition. It should be noted that the dominant energy 
condition holds for F being the power function of / or J, whereas it is not the case when 
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F is given as a trigonometric function of its arguments. Note that in presence of A-term 
the role of other parameters such as order of nonlinearity n, perfect fluid parameter C, and 
spinor mass in the evolution process are rather local, while the global process are totally 
determined by the A-term, e.g., for a positive A we have always oscillatory mode, while 
for a negative A solutions are generally inflation-like though for some special choices of 
problem parameters the oscillatory mode of evolution can be attained. 
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